Periodic Coherence Peak Height Modulations in Superconducting Bi 2 Sr 2 CaCu208+,5 
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In this paper we analyze, using scanning tunneling spectroscopy (STS), the local density of elec- 
tronic states (LDOS) in nearly optimally doped Bi2Sr2CaCu20g+i in zero field. We see both 
dispersive and non-dispersive spatial LDOS modulations as a function of energy in our samples. 
Moreover, a spatial map of the superconducting coherence peak heights shows the same structure 
as the low energy LDOS. This suggests that these non-dispersive LDOS modulations originate from 
an underlying charge-density modulation which interacts with superconductivity. 

PACS numbers: 74.72.Hs, 74.50.+r, 74.25.-q 



I. INTRODUCTION 

The Scanning Tunneling Microscope (STM) has been 
an important tool in the study of high-temperature su- 
perconductors since their discovery. Initially, a variety 
of gap sizes and structures were found and introduced 
much controversy into the subject. Later, a more coher- 
ent consensus among different groups emerged regarding 
the surface properties of these high-T c materials. To give 
a few examples, STM studies revealed the nature of the 
superstructure in BSCCOi, the d-wave nature of the gap 
and its size^, the effect of local impurities, the emergence 
of zero-bias anomalies 3 *^, and the electronic structure of 
the core of vortices^. More recent measurements sug- 
gest that superconductivity may not be homogeneous in 
high-T c superconductors. In particular, STM measure- 
ments have found spatial variations of the gap size in 
YBCO& and BSCCO 9 i 10 i n . 

While gap inhomogeneities have been found to dom- 
inate the electronic structure at large measured bias, 
more ordered structures underlying the d-wave-like tun- 
neling spectra have been found at lower energies. A 
current topic of great interest in high-T c supercon- 
ductors is the presence of spatial modulations of the 
charge and spin densities. Theoreticali2ii 3 riiiii5iiSiil 
and experimental 18 i 19 i 20 i 21 i 22 i 23 i 24 i 25 evidence has been 
mounting in support of the possibility that their ground 
state exhibits spin and charge density waves (SDW and 
CDW), which may be primarily one-dimensional (i.e., 
stripes) , or two-dimensional2& with a characteristic wave 
vector in the Cu-0 bond direction of q^-o = 0.25(27r/ao). 
In STM measurements such a modulation was first seen 
by Hoffman et alSL in a magnetic field as a 2D check- 
board pattern of LDOS, aligned with the Cu-0 bonds, 
around vortex cores in slightly over-doped BSCCO with 
Ni impurities. The reported modulations showed a 
checkerboard ordering vector of q^-o = 0.23(27r/a ) ex- 
tending to large distances when measured at bias en- 
ergy ~ 7 meV. Howald et alm& shortly afterwards re- 



ported this same effect in zero field on similarly doped 
BSCCO crystals without intentional substitution impuri- 
ties. The observed modulation with ordering wave vector 
q-n-o = [0.25±0.03](27r/ao) was found at all energies, ex- 
hibiting features characteristic of a two-dimensional sys- 
tem of line objects. Moreover, Howald et al. showed that 
the LDOS modulation manifests itself, for both positive 
and negative bias, as a shift of states from above to below 
the superconducting gap. The fact that a single energy 
scale (i.e. the gap) appears for both superconductivity 
and these modulations suggests that these two effects are 
closely related. 

In subsequent studies at zero field, Hoffman et alM> and 
McElroy et almS. measured the dispersion of the strongest 
Fourier peak along the it — (i.e. Cu-O) direction. They 
asserted that it was consistent with what is expected from 
quasiparticle scattering interference 31 , in which a peak in 
the Fourier LDOS equals the momentum transfer wave 
vector of the incident and scattered waves. In general, 
their data showed good agreement with photoemission 
results (i.e. band structure results 3 ^) at large bias, but 
did not continue to disperse below ~ 15 mV. 

To account for all the available data in the full en- 
ergy range, Vojta- 3 ., Podolsky et al~ and Kivelson et 
alm^ have shown that both charge order and quasiparti- 
cle scattering effects can occur in the presence of pinned 
fluctuating stripes. In particular, Podolsky et al. , using 
explicit calculations, showed that for a system with incip- 
ient charge order the dispersion for this particular effect 
is very weak and less than expected by band-structure. 
At low energy it converges to the ordering vector, rather 
than the vector corresponding to the nodal separation 
on the Fermi surface. A similar conclusion was reached 
by Kivelson et al. who emphasized that at higher en- 
ergies, for a relatively clean material, quasiparticle scat- 
tering interference will show a strong signal that over- 
whelms the weak charge modulation. At low energies, 
the minimum energy required to overcome the finite su- 
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perconducting gap means that only quasiparticles near 
the nodes participate and thus should give a wave vec- 
tor for the interference which is larger than what was 
measured. The discrepancy was taken as evidence that 
another phenomenon dominates at low energies. 

In this paper, we present new data and analysis on this 
phenomenon, measured on near optimally doped sam- 
ples. Using the same apparatus with the same exper- 
imental conditions as before, we show data for two of 
our samples (from the same growth run): one which has 
a (fixed) two-dimensional ordering wave vector g T _o = 
[0.25 ± 0.03](2-7r/ao) that dominates the dispersive sig- 
nal over much of the energy ranged, and one which 
has a stronger dispersive signal that begins to dominate 
from lower energies on up, yet still has a majority of 
spectral weight for the lowest energies at wave vector 
q^-Q = [0.22 ± 0.03] (2n/a ). We further show that for all 
the samples studied, the density of states at the gap (i.e. 
dl/dV at V = A) exhibits modulations with a wave vec- 
tor similar to the low-energy pattern. Finally, we present 
evidence for the interplay of these modulations and su- 
perconductivity by showing the strong supression of the 
large gap coherence peaks. 



II. EXPERIMENT 

We performed measurements using a home-made cryo- 
genic STM. The STM measures differential conductance 
G = dl/dV which is proportional to the LDOS. The sam- 
ples are near optimally doped (slightly over-doped) single 
crystal Bi 2 Sr 2 CaCu20 8+< 5 (T c - 86-87 K)^ 6 . grown by a 
floating-zone method. They are cleaved at room temper- 
ature in an ultra high vacuum of better than 1 x 10~ 9 
Torr, revealing an atomically flat surface between the 
BiO planes. Then they are quickly lowered to the cryo- 
genic section at a temperature of 6-8 K, where cryopump- 
ing yields orders of magnitude better vacuum. Typical 
data were taken with a sample bias of —200 mV and a 
set point current of —100 pA, which establishes the rel- 
atively arbitrary normalization for the LDOS. We also 
performed measurements with a sample bias of +65 mV 
and a setpoint current of +25 pA. 

At each point on the surface, a spectrum (dl/dV vs. 
sample voltage V) was taken. The bias modulation for 
the spectra is 1 mV rms . This, in addition to the time 
constant of the lock-in amplifier used to record dl/dV, 
yields a total blurring of the spectra of ~ 3 mV. Although 
there are discrepancies when trying to fit a spectrum with 
a d-wave function, using the voltage, i.e. position, of the 
coherence peak maximum for the gap value A yields a 
reasonable fit, and is the method we will use throughout 
this paper. The coherence peak-heights map is made 
by evaluating dl/dV at V = A (at that location). All 
maps have been constructed for positive sample voltage, 
as this yields better signal to noise. This observation is 
most likely related to the asymmetry of the conductance 
spectra which are common to all BSCCO STM studies^!. 



III. RESULTS AND ANALYSIS 
A. Spatial Variation 

Spatial variations of the superconducting gap on the 
surface of BSCCO have been reported by several groups 
and their existance is now an established fadpi2*ii*2&. 
Typically the size of the gap varies on a length scale 
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FIG. 1: a) Gap size A over a 140 A x 140 A area, b) Coher- 
ence peak height (divided by the average conductance) with 
inverted colorscale. Scan performed at +65 mV bias. 

of roug hly 30 A as was described in detail elsewhere by 
Howald et aJiS. The gaps vary from small gap regions of 
~ 30 mV to large gap regions of ~ 60 mV. The regions of 
large gap usually develop from average gap background 
with the largest gap at their center—. Fig. shows the 
size of the gap mapped over a typical surface. Patches 
about 30 A across of varying gap size are clearly visible. 
The distribution of gap sizes in this view is depicted in 
Fig. [2Ji, where the average gap is A = 45.4 mV. The 
smallest scale features reflect some variation with atomic 
resolution. In addition, the partial, nearly vertical lines 
show that there is some correlation between superstruc- 
ture and the gap (Fig.[IJi). 

An important feature of the gap-size distribution is 
that the height of the coherence peaks is varying as well. 
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FIG. 2: a) Histogram of gap size distribution for the data in 
Fig. b) Density plot of coherence peak height versus gap 



There is some correlation between these two effects: those 
spectra with the smallest gaps typically have taller coher- 
ence peaks. This anti-correlation is depicted in Fig. |2t>, 
which shows the variation of the gap size and the coher- 
ence peak height for the area in Fig. ^ The colorscale 
of Fig. is inverted in order to demonstrate the anti- 
correlation shown in Fig. 03 (e.g. regions with large gaps 
tend to have shorter coherence peaks). 

While there is a correlation between gap size and co- 
herence peak height, it is not a simple one-to-one re- 
lationship. Thus the two maps can differ significantly. 
Aside from some high frequency noise, the spatial vari- 
ation in the peak height (Fig. is smoother, with no 
abrupt transitions between the regions. Next, we note 
the presence of an atomic corrugation which is stronger 
in one direction. This is most likely due to a convolution 
with the tip shape, as this effect shows up in the topo- 
graphic image as well. More importantly, it is clear that 
the peak-height shows an ordered structure (see e.g. the 
lower-right corner of the hgure) in which the modulation 
amplitude can be as great as ~30% of the mean peak 
height in some areas. 

By inspection of spectra from —200 mV to +200 mV, 
we found that normalization by the setpoint current at 
+65 mV makes the peak height map least sensitive to 
contributions from the superstructure as well as most of 
the low fc-vector structure (which is due to gap size inho- 
mogeneities) . Additionally, choosing a positive normal- 
ization voltage makes the spectra less dependent on the 
normal state background which is typically stronger on 
the negative bias side of the gap^I. Thus taking the scan 
at +65 mV bias largely removes these features, as com- 
pared to a peak height map taken with a more common 
setpoint voltage of —200 mV. This simple procedure al- 
lows the peak height structure to be seen in the real space 
image. A map of the current at +65 mV, for a spectro- 
scopic scan taken at —200 mV, shows little to no spectral 
weight above the noise near q = 0.25(27r/ao), and thus 
this procedure would not create the modulations we see. 



B. LDOS Modulations, Dispersion, and G(A) 

To look for LDOS modulations, one typically looks at 
the differential conductance G = dl/dV as a function 
of voltage (V) at each point on the sample. In most 
cases modulations are more visible in Fourier space where 
the length scales of various features are better separated. 
Fig. |3 shows a Fourier transform of the differential con- 
ductance of the area shown in Fig. ^ for two different 
energies (10 mV and 29 mV) as well as for the gap, A 
(note that the majority of gaps are 30 mV < A < 60 mV 
for that region). Circles are placed at (27r/ao)(±0.25,0) 
and (27r/ao)(0,±0.25) as reference points. 

As in previous results^ 8 (Fig. , for low energies we 
sec periodic density of states modulations at a periodicity 
close to 4eto. The periodicity (as shown by the Fourier- 
analysis) is [0.22 ± 0.03](27r/a ) for the data in Fig. 
At low energies (< 15 mV), these modulations dominate 
(Fig. |3^) , but as one goes up in energy, the strength of 
the overall signal increases and moves to longer wave- 
lengths in a very similar way to the results of McElroy 
et alm&. This behavior is clearly seen in Fig. |3Jd where 
we show a Fourier transform for a sample bias of 29 mV. 
(The color scale on the Fourier maps has been adjusted to 
keep this signal within view.) At higher energies, above 
~ 35 mV, the signal is lost in the noise due to the appear- 
ance of coherence peaks that vary with position on the 
sample. Although our experimental fc-resolution limits 
precise quantitative statements we can make about any 
dispersion, it is clear that overall, more spectral weight 
is appearing at lower fc-vectors as the energy increases. 

To show this dispersive effect, we take line scans along 
the 7T — direction in Fourier space (Figs. EJi and^Jl). 
At each point along the line, we weight the neighbor- 
ing Fourier points with a Gaussian filter of FWHM 1.6 
pixels (which is not necessarily centered directly over a 
pixel). The values are then squared, summed, and a 
square root is taken. Finally, this value is normalized by 
the proximity of the line to the various Fourier points. 
The values reported in the line scans are the total mod- 
ulation amplitude in that region of Fourier space. Due 
to pixelation effects and the size of our scans, the uncer- 
tainty in the peak positions, as well as our resolution, is 
Aq = 0.03(2-7r/a ). However, the normalization proce- 
dure and large filter width ensure that all of the features 
(i.e. peaks) seen are distinct and not due to pixelation 
artifacts (e.g. a broad peak will not be split into two by 
passing the line scan between pixels). 

We observe that at all the energies shown a (sometimes 
weak) signal always exists near the four-period wavevec- 
tor, in agreement with Howald et al. 28 . Additionally, 
there exists a signal at a slightly lower /c-vector. This is 
most likely the dispersive quasiparticle scattering inter- 
ference peak since it becomes suppressed upon an inte- 
gration over energy, as shown for example in the above 
work2&, Fig. 9. If one follows the "main peak" as sug- 
gested by Hoffman et al. (i.e. using a single peak to 
fit both features) then dispersion is inferred at energies 
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FIG. 3: a) FFT of LDOS (dl/dV) at 10 mV. White arrow FIG. 4: Previously published data, a) FFT of LDOS at 11 

indicates direction of line scan, b) FFT of LDOS at 29 mV c) mV. b) FFT of LDOS at 28 mV c) FFT of LDOS taken at 

FFT of LDOS taken at the coherence peak maximum d) Dis- the coherence peak maximum with 65 mV normalization (see 

persion relation of the charge modulation periodicity. Black Section C.) d) Dispersion relation of the charge modulation 

line is coherence peak maxima. All lines shifted for clarity. periodicity. Black line is normalized coherence peak maxima. 
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above ~15 mV. 

If we consider the modulations in the peak heights 
(Fig. njj) in Fourier space (Figs. Bt,d), we find that it 
is very similar to the low energy modulation (Fig. yk>). 
The peak-heights map has a periodic structure close 
to four lattice spacings (more precisely q-^-o = [0.22 ± 
0.03](27r/ao)) even though the entire contribution comes 
from energies above ~ 30 mV (see Fig. By looking 
at the LDOS based on gap size, we are sampling from 
a range of energies whose spectra contribute in such a 
way as to give a modulation at a higher q than would be 
predicted by quasiparticle scattering for any of these en- 
ergies alone. The peak (or shoulder in the case of Fig. 0Ji) 
at slightly lower q is likely to be from these quasiparticles 
at higher energies. 

Finally, we note that to maximize the amplitude of the 
peaks, our line scan in Fig. [JJi is rotated by ~4 degrees 
with respect to the atomic positions. This is not sur- 
prising considering that our scan size is only over a few 
correlation lengths, and thus defects can cause an overall 
rotation in the modulations. 



C. Interaction with Superconductivity 

Some other interesting effects can also be seen when 
one compares the maps of the coherence peak heights and 
the gap size. Here we show previously published data2& 
that has been taken over a larger area (Fig. O ■ A novel 
procedure we are using here is to normalize (i.e. divide) 
the individual spectra by the current at +65 mV. (More 
precisely, we divide by the average conductance from to 
65 mV.) This is similar to scanning over the area with this 
setpoint voltage of 65 mV (as in our previous figures) . We 
note that this procedure is only used to "clean up" the 
image so that certain features can be clearly seen in real 
space (note the lack of low frequency noise in Fig. 0k). 
As mentioned earlier, this procedure does not add any 
artificial modulations. 

By comparing Figs.^ andJSJj (or Figs.^ and^>), one 
can see that the amplitude of the coherence peak DOS 
modulations is larger in the regions of large gap. In con- 
trast, regions of small gap show modulations of reduced 
amplitude. Since there are only a few modulation crests 
and troughs within a particular region of large or small 
gap, this effect is difficult to quantify, although it can 
most easily be seen by following the regions of largest 
gap. We note that the large amplitude modulations are 
not simply due to the normalization procedure (regions 
of large gap are renormalized to have higher LDOS due to 
the incomplete integration of the coherence peaks), since 
they are actually a larger fraction of the mean coherence 
peak height. 

Fig. |SJ: shows the real space LDOS at 8 mV, af- 
ter Fourier filtering near the (27r/ao)(±0.25,0) and 
(27r/a )(0,±0.25) points in reciprocal space. We chose an 
energy where the signal from the underlying order domi- 
nates over the quasiparticle scattering interference signal. 
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FIG. 5: a) Gap size distribution over a 160 A xl60 A area 
b) Coherence peak heights, normalized c) LDOS at 8 mV, 
un-normalized, Fourier filtered. Equivalent points on Figs, b 
and c are marked by points A, B, C. 



The filter is shaped like a circle, with a radius of 4 pix- 
els in which no Fourier weight is suppressed. The edges 
then taper off like a Gaussian with a width of one pixel. 
This filter encompasses the vast majority of the spectral 
weight in the region and thus does not favor one partic- 
ular wave vector. The filtered image shows a dominant 
four-period modulation that is almost checkerboard like, 
but with dislocations in the form of extra half-rows. By 
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following the modulations, one can see a general corre- 
spondence between the peaks of Fig. [5Jd and the troughs 
of Fig. Et, and vice versa. The points marked A, B, 
and C are examples of this out-of-phase relation. (The 
features may not match up exactly due to noise or the 
slight contribution from quasiparticle scattering.) This 
suggests that where the low energy modulations have an 
increased LDOS, the coherence peaks are suppressed and 
vice versa. In particular, this suppression is stronger for 
the large gap coherence peaks. 

IV. DISCUSSION 

It has been argued before that the competition be- 
tween kinetic energy and Coulomb repulsion may lead 
to various forms of charge and spin ordered states. In 
particular, "stripes" have been predicted to occur in 
doped antiferromagnet9i2i±2i22iffi. The discovery^i* 4 ^ 3 ^ 
of stripe order in L^-xSr^NiC^+i and soon after in 
Lai.6-xNdo.4Sra;Cu04^ added considerable credibility 
to the suggestion that stripe states form an important 
bridge between the Mott insulator and the more metallic 
state which becomes Fermi-liquid-like at heavy doping. 
The theory is that at low doping, static stripes char- 
acterize a true broken symmetry state. However, upon 
increased doping the periodicity of the stripes decreases 
and coupling occurs. This is roughly the region where su- 
perconducting order wins over stripe order. Stripes will 
now exist in the superconducting phase in a dynamical 
sense, i.e. stripe order is fluctuating. Therefore no broken 
symmetry is expected and thus to preserve the point sym- 
metry of the underlying crystal, a fluctuating checker- 
board is expected. STM, however, is a static probe and 
thus cannot detect any structure associated with fluctu- 
ating order unless something pins itii^SiLiiiSL22i2Si^. 
Indeed, the gap size inhomogeneities discussed above and 
other forms of chemical disordei^ are a natural source 
for pinning and can make stripes or checkerboards visible 
to STM in the form of LDOS modulations. 

Four-fold symmetrical order may also be a consequence 
of strong interactions on the square lattice. When re- 
duced to the low energy Plaquette Boson Fermion Model, 
the system shows a checkerboard structure due to the 
tendency of this model to locally prefer a four boson (an 
antifcrromagnon triplet and a d-wave hole pair) state (i.e. 
doping of 1/8) 26 47 . In cither of the above theories, this 
underlying order has a non- or weakly dispersive nature. 
In the case of "striped" structures, weak dispersion is ex- 
pected due to the finite size of the stripe domains and 
the interaction with the itinerant quasiparticlesi^ 4 ^ 

While the above two examples result in an (almost) 
"fixed-q" order (i.e. a true frozen charge density visi- 
ble in the LDOS at all energies at an almost constant 
wave vector) , the existence of quasiparticles (in the broad 
sense 35 ) in the presence of weak disorder may also add 
quasiparticle scattering interference effect a 31 i 48 . In this 
case, quasiparticles of given energy scatter off an im- 
purity. The resulting interference between the origi- 



nal and scattered waves leads to variations of the local 
density of states at wave vectors q = k k', where 
k and k' are the wave vectors of states with energy 
E = e(k) = e(k'), as determined by the band structure, 
e(k). Judging from measured band structure on BSCCO 
via photoemissionS, quasiparticle scattering interference 
effects should be strongly dispersive, as is indeed seen in 
STM experiments for energies greater than ~ 15 mV. 

However, analysis of all recent experiments 3 ^ indicates 
that both evidence for a fixed-g oscillation^ and quasi- 
particle scattering interference22iSI have been found ex- 
perimentally. Taking this point of view, it is clear that 
if there is an underlying order coexisting with the quasi- 
particle interference structure, one needs a way to sep- 
arate out these effects. The main problem here is that 
the large contribution of the gap inhomogeneities and the 
strong dispersion of the quasiparticle scattering interfer- 
ence cover up this underlying order. Following the ideas 
of Kivelson et ai, Howald et al. showed that integration 
of the Fourier space LDOS over a wide range of ener- 
gies reduces the influence of any random or dispersing 
features while at the same time enhances features that 
do not disperse. Another possibility is to look for the 
interaction of these features with another order parame- 
ter, such as superconductivity. We claim that the peri- 
odic structure observed in the LDOS at the gap, G(A) 
in Fig. \I]p, is exactly this effect. 

What we see is a return of features from the low en- 
ergy region (where saturation is observed in the plot of 
dispersion), suggesting that these features indeed exist 
and are related to superconductivity. In addition, the 
large amplitude of these coherence peak modulations and 
the distinct peak near q — 0.22(27r/ao), as shown in the 
line scan of the Fourier transform of the coherence peak 
heights (Figs. Oil and^Jl), is suggestive of some kind of co- 
herent contribution from a non-dispersive feature. (Dis- 
persive contributions would cause a peak in the line scan 
to be broadened or suppressed, e.g. the peak/shoulder 
at lower k- vector.) 

Our Fourier analysis from low energies up to the small- 
est gap sizes (where the noise from inhomogeneities over- 
whelms our signal) supports the picture presented ear- 
lier of a non- or weakly dispersive feature near q = 
0.22(27r/ao) in addition to a dispersive feature at a lower 
A;-vector. We find that the large amplitude of this lower 
k-vectov feature swamps out the non-dispersive feature 
at higher energies, but becomes weaker as one goes down 
in energy, until at approximately 15 mV only the non- 
dispersive signal remains, thus explaining the saturation. 
The additional line scan for the Fourier transform of co- 
herence peak heights can be seen as a way to remove the 
effects of gap size inhomogeneities to reveal that a struc- 
ture at q « 0.22(27r/ao) still exists at higher energies. 
It makes a similar point as the spatial maps of the co- 
herence peaks, namely, that by selectively sampling from 
the (higher) energies related to superconductivity the 
low energy features reappear. 

Finally, we note the correspondence of these (coherence 
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peak height) modulations to the superconductivity as de- 
termined by the size of the gap. First we observe that the 
modulation structure that resides in the regions of larger 
gaps is the most pronounced. These regions of large gap 
with low coherence peaks generally resemble slightly un- 
derdoped samples 10 . On the other hand, the modulation 
is suppressed and the coherence peak heights are more 
uniform in regions of small gap and tall coherence peaks. 
Gaps in these regions are more similar to gaps found in 
overdoped samples^ 9 -; this is consistent with the notion 
that beyond optimal doping a more homogeneous charge 
density exists as the system tends more towards a Fermi 
liquid state. 

Such an observation does not necessarily point out 
a competition between charge-density modulation and 
superconductivity, but rather reinforces the idea that 
the two effects coexist at and below optimal dop- 
ing. One possible interpretation is that the fluctuat- 
ing stripe/checkerboard phase exists in all the regions 
below optimal doping, and as one moves further into 
overdoping (i.e. into regions of small gap) the modu- 
lations become diminished. Our observations therefore 
complement those of Vershinin et alm& who found simi- 
lar patterns in the pseudogap regime of slightly under- 
doped Bi2Sr2CaCu208+5. As noted by Kivelson et alj&, 
the effect of quasiparticle scattering interference should 
disappear at temperatures above T c , revealing the under- 
lying order. For our measurements at low temperature, 
in regions of very large gap with weak coherence peaks 
(similar to the pseudogap), charge ordering is indeed visi- 
ble. The two results therefore suggest that in the absence 
or suppression of superconductivity, charge-ordering may 
be the preferred phase. 

However, these "pseudogap-like" regions are only 
about two superconducting coherence-lengths large, 
while the modulations (as seen in Fig. (5J;) persist over 
approximately seven periods^SS, suggesting that these 
modulations are more apparent in, but not exclusive to, 
these regions. Moreover, both the modulation and su- 
perconductivity seem to coexist at low temperatures for 
all gap sizes except that in the small gap regions the 



peak-height modulation amplitude is suppressed as dis- 
cussed above. However, since the pseudogap is likely 
not a true phase transition, and in terms of both dop- 
ing and temperature the system may be relatively far 
from a charge ordering critical point, ordering may be 
tenuous. Upon lowering the temperature, the system un- 
dergoes a true phase transition into the superconducting 
phase where one expects the fixed-wavelength modula- 
tion phenomenon to be strongly pinned in the presence 
of disorder. We see this effect as a non-dispersive signal 
which additionally manifests itself in the coherence peak 
heights of the tunneling spectra. 



V. CONCLUSION 

The observation of spatial modulations in a quan- 
tity that is related to the pair amplitude, i.e. the co- 
herence peak heights, points to the intimate relation 
between superconductivity and the charge modulation. 
This is strong evidence that the underlying modulation 
at q w 0.22(27r/ao) which exists down to low energies, 
yet reappears at gap energies, is a separate effect from 
the quasiparticle scattering interference signal. Further 
analysis of the energy dependence suggests that the "sat- 
uration" seen in plots of dispersion can be explained by 
the dominance of this non-dispersive signal at low en- 
ergies. Finally, the observation that these modulations 
suppress the coherence peaks most strongly in regions of 
large gap is consistent with a phase diagram that has an 
ordered phase towards under-doping and tends towards 
a more homogeneous charge distribution with increased 
doping. 
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